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1. Introduction 

Quark colour confinement in hadron physics remains unexplained after more than 30 years of 
intense study (for a recent overview see Ref. [1]). In lattice gauge theory two avenues of research 
have been most commonly adopted: confinement by Z{N) centre vortices and confinement by 
means of abeUan monopoles (for a critical discussion of both see Ref. [2]). Gauge fields can be 
first fixed to a particular gauge, such as Maximal Abelian Gauge (MAG) or Maximal Centre Gauge 
(MCG), monopoles and centre vortices are then defined by the projection of these gauge-fixed fields 
to U{1)^^^ or Z{N) gauge fields, respectively. Much of the progress to date has occurred in SU (2) 
using MAG and MCG, with original findings reproducing about 90% [3] and about 100% [4] of the 
non-abelian string tension. Removing monopole [5, 6, 7] or centre vortex [7, 8, 9, 10, 1 1] degrees 
of freedom from SU (2) lattice gauge fields appears to leave topologically trivial, non-confining 
gauge fields that preserve chiral symmetry. 

Work in SU{3) has not progressed to this level. While initial investigations were hopeful 
[12, 13], subsequent results for MAG and MCG were not so encouraging [14, 15], both failing to 
reproduce the full non-abelian string tension. Earlier studies in SU (2) using MCG showed that 
the centre-projected configurations recovered the full string-tension, however further study into 
the ambiguities of the gauge-fixing procedure showed that this result is plagued by Gribov copy 
effects [16, 17]: methods which give higher values of the gauge fixing functional produce smaller 
values for the vortex induced string tension. We point out that when the Laplacian Centre Gauge of 
Ref. [18] (which is free of Gribov ambiguities on the lattice) is used as the fixing method the full 
SU{3) (and SU (2)) string tension is recovered for the centre-projected gauge fields. However, the 
interpretation of this vortex matter is cumbersome in the continuum limit [15, 19]. 

In this paper we focus on the Gribov problem of the SU (3) centre vortex picture of confinement 
using the MCG gauge-fixing method. In the centre- vortex picture the gauge fields are considered to 
be decomposed into a long-range field carrying all the confining fluctuations and a short-range 
field containing non-confining perturbations as well as other short range effects 

Uf,{x)=Z^{x)V^{x). 

Here Z^(x) is the centre element which is closest, on the SU (3) group manifold, to U^{x). 

1.1 Smearing as a Preconditioner 

Since the centre elements, in the centre vortex picture of confinement, correspond to the long- 
range physics, we employ the use of smearing, which smoothes out the short-range fluctuations, to 
construct a preconditioning gauge transformation for each gauge field prior to gauge-fixing [20]. 

Firstly, we smear the gauge field using any smearing algorithm (stout-Unk smearing [21] was 
used in the data shown here but other smearing algorithms are currently under investigation). We 
then fix the smeared field using the MCG gauge-fixing method. 

At each iteration we keep track of the total gauge transformation that has been applied to the 
smeared gauge field. Once the algorithm has converged we use the stored total transformation 
as a preconditioning gauge transformation for the unsmeared gauge field. We emphasise that the 
(unsmeared) preconditioned gauge field remains on the same gauge orbit since the preconditioning 
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is merely a (specific) gauge transformation on the original links. Gauge-fixing the preconditioned 
field simply gives us a Gribov copy of the result from gauge-fixing the original gauge field. 



2. Identifying Vortex Matter 

We employ the MCG gauge-fixing algorithm as outlined in Ref. [15]. The gauge condition we 
chose to maximise (with respect to the centre elements Z^(x)) in this algorithm is 



Vum = l£ [itr[/,"(x)] [\tvUf^{x)]\ 



where Ni is the number of links on the lattice and is the gauge-transformed field. 

After fixing the gauge, each link should be close to a centre element of SU{3), Z'" = e^^'" , 
(j)'" = with m e { — 1,0, 1}. Since, for every Unk, 

-trUf^ (x) = M,. . e'^-" and (/)^ . = tan" > ^ — 

3 ^ ^ Re(tri7^(x)) 

then ^x,^i should be close to some 0"', by construction of the gauge-fixing condition. We then 
perform the centre projection by mapping 

SU{N)^Zn: U^{x)^Z^{x) with Z^{x) = e'^^-i^ , 

with the appropriate choice of 0"'^, m G { — 1,0, 1}. 

To reveal the vortex matter we simply take a product of links around an elementary plaquette. 
We say a vortex pierces the plaquette if this product is a non-trivial centre element and the plaquette 
is then a P-vortex. In a toy Yang-Mills model we can remove these P-vortices by hand from the 
configuration using U'^{x) = zJi{x)U^{x). 



3. Results 



Calculations are performed using 100 quenched configurations with the Luscher-Weisz pla- 
quette plus rectangle gauge action [22] on a 16^ x 32 lattice with j3 = 4.6. Similar results are 
being found on 200 20^ x 40 lattices but since these results are currently incomplete they are not 
presented here. 

Stout-link smearing with a smearing parameter of 0.1 is used to construct the preconditioning 
transformation with the number of sweeps ranging from to 12 in steps of 4 sweeps. Here, each 
preconditioning was conducted independently although it would be possible to use the precon- 
ditioning from an initial level to precondition subsequent levels of smearing, thereby decreasing 
computation time when seeking comparative results. 

If we look at the number of P-vortices located in each configuration compared with the number 
located on the next level of preconditioning (Table 1), we can see that the first level of smearing 
decreases the number of vortices identified in all configurations. This behaviour continues right 
through all levels of preconditioning although the size of the effect drops to 79% of configurations 
by the last preconditioning. The magnitude of the reduction in number of P-vortices is as high as 
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Comparison 


Decreased 


Higher 


Sweeps 


Vortices 


Maxima 


0^4 


100% 


100% 


4^8 


85% 


79% 


8^ 12 


79% 


64% 



Table 1: The percentage of configurations 
that have a lower number of P-vortices 
and the percentage of configurations that 
achieve a higher gauge condition maximum 
when comparing two different values of the 
preconditioning. 



Sweeps 


Total Iter. Blocks 


Smear Max 


Max 





60 ±16 


**** 


0.7360(11) 


4 


89 ±15 


0.9153(16) 


0.7409(9) 


8 


97 ±20 


0.9375(21) 


0.7417(11) 


12 


95 ±17 


0.9461(17) 


0.7421(9) 



Table 2: For each of the sweeps used in the precondi- 
tioning: the average total (smeared gauge field fixing 
plus preconditioned gauge field fixing) number of iter- 
ation blocks (of 50) used, the average smeared gauge 
condition maximum reached and the average precon- 
ditioned gauge condition maximum reached. 



66 ± 16% when using the first level of preconditioning with the effect dropping to as low about 
10% for the transition between the 8 and 12 sweep smearing levels. The gauge condition maxi- 
mum achieved also increases in all configurations for the first level of preconditioning, with the 
magnitude of this effect dropping to 64% by the final preconditioning. 

We also note from Table 2 that the total number of blocks of 50 iterations of the gauge-fixing 
algorithm increases slightly when using the preconditioning. Typically, |'s of the iterations are 
spent fixing the smeared field and j fixing the preconditioned field. The gauge-fixing maximum 
achieved monotonically increases however, both for the smeared fields and for the preconditioned 
field (though not to the same degree) with increasing levels of smearing preconditioning. It should 
be noted that, regardless of the preconditioning level, the centre phases of the links of the fields 
always remain evenly distributed across the three possible values, reflecting the fact that the reali- 
sation of centre symmetry remains unaffected. 

3.1 The Static Quark Anti-quark Potential 

Computing the static quark anti-quark potential as a function of their separation is a two step 
process. Wilson loops W{R,T) of extension of RxT have the large T behaviour 

{W{R,T)) ocexp{-V {r)aT}, r:=Ra, 

where a is the lattice spacing. The method for extracting the potential is identical to that of two- 
point functions. To obtain the static quark anti-quark potential as a function of the quark separation 
we simply repeat this process for a range of values of the separation R. As can be seen in Figure 1, 
by using off-axis spatial paths for the Wilson loops, we can obtain non-integer values of R. 

In Figure 2 the effective potential V{r) is plotted against T (by taking the log of the ratio 
of the value of the Wilson loops at two adjacent time slices) for two Gribov copies of each of the 
vortex-only and vortex-removed configurations and comparing these to the results from the original 
configuration. On the left is the unpreconditioned Gribov copy and on the right is a Gribov copy 
that uses 4 sweeps of smearing. We are seeking to find a plateau in the potential in each case. 
We note that the quality of the fit, as defined by the X'^ values, depends greatly on the fit range 
chosen and one must first account for the systematic drift in the data before selecting a sensible fit 
region. In particular, for the vortex-removed configurations, the potential does appear to become 
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degenerate at larger values of the separation but choosing a suitable fit is challenging. We also 
note that the scale of the potential is greatly reduced for the vortex-only configurations while the 
reliability of the result increases dramatically, allowing for fits at very large T values. 
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Figure 1: The static quark anti -quark potential plots for each of preconditioning smearing sweeps used. 
Each plot contains data for the full, vortex-removed and vortex-only configurations. 



In Figure 1 the potential V{r) is plotted against the separation r for our original configuration 
and 3 Gribov-copy ensembles. We can see that the results are consistent with the loss of confine- 
ment for the vortex-removed configurations. The preconditioned configurations better reproduce 
the full potential at small distances. We also note that the string tension of the vortex-only results 
drops dramatically when we apply our preconditioning, going from about 65% of the full string 
tension to as low as about 26%. 



4. Conclusions 

We have shown that using smearing as a preconditioner to the MCG gauge fixing algorithm in- 
creases the maximum of the gauge-fixing condition. We find that these higher maxima correspond 
to a lower number of measured of P- vortices. When looking at the static quark anti-quark potential 
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for the vortex-only and vortex-removed configurations we find results consistent with loss of con- 
finement for all the vortex-removed configurations. We also find that the preconditioning decreases 
the measured string tension of the vortex-only results, reducing them from 65% to as low as 26% 
of the full string tension. 

Similar to what has been observed in SU{2) [16, 17], there appears to be a significant anti- 
correlation between the value achieved in the gauge fixing functional and the percentage string 
tension reproduced by centre vortices. Although the fundamental modular region of MCG would 
be an ideal candidate for a unique definition of vortex texture, it seems that the vortex matter arising 
from the first Gribov region as a whole has a bigger phenomenological relevance. Further work is 
continuing on larger lattice volumes. Results to date are consistent with findings here. 
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Figure 2: The effective potential plots for each of the untouched, vortex-removed and vortex-only configu- 
rations. On the left are the results from the unpreconditioned gauge-fixed fields and on the right the results 
from preconditioned gauge-fixed fields with 4 sweeps of smearing used in generating the preconditioning 
gauge transformation 



